Increasing thermoelectric efficiency towards the Carnot limit 
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We study the problem of thermoelectricity and propose a simple microscopic mechanism for the increase 
of thermoelectric efficiency. We consider the cross transport of particles and energy in open classical ergodic 
billiards. We show that, in the linear response regime, where we find exact expressions for all transport co- 
efficients, the thermoelectric efficiency of ideal ergodic gases can approach Carnot efficiency for sufficiently 
complex charge carrier molecules. Our results are clearly demonstrated with a simple numerical simulation of 
a Lorentz gas of particles with internal rotational degrees of freedom. 
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Although thermoelectricity was discovered almost 200 
years ago, a strong interest of the scientific community arose 
only in the 1950's when Abram loffe discovered that doped 
semiconductors exhibit relatively large thermoelectric effect. 
This initiated an intense research activity in semiconductors 
physics which was not motivated by microelectronics but by 
the loffe suggestion that home refrigerators could be built with 
semiconductors [1, 2]. As a result of these efforts the thermo- 
electric material Bi2Te3 was developed for commercial pur- 
poses. However this activity lasted only few years until the 
mid 1960's since, in spite of all efforts and consideration of 
all type of semiconductors, it turned out that thermoelectric 
refrigerators have still poor efficiency as compared to com- 
pressor based refrigerators. Nowadays Peltier refrigerators 
are mainly used in situations in which reliability and quiet 
operation, and not the cost and conversion efficiency, is the 
main concern, like equipments in medical applications, space 
probes etc. In the last decade there has been an increasing 
pressure to find better thermoelectric materials with higher 
efficiency. The reason is the strong environmental concern 
about chlorofluorocarbons used in most compressor-based re- 
frigerators. Also the possibility to generate electric power 
from waste heat using thermoelectric effect is becoming more 
and more interesting [1-3]. 

The suitability of a thermoelectric material for energy con- 
version or electronic refrigeration is evaluated by the thermo- 
electric figure of merit Z, 



(1) 



where a is the coefficient of electric conductivity, S is the 
Seebeck coefficient and k is the thermal conductivity. The 
Seebeck coefficient S, also called thermopower, is a measure 
of the magnitude of an induced thermoelectric voltage in re- 
sponse to a temperature difference across the material. 

For a given material, and a pair of temperatures Th and Tq 
of hot and cold thermal baths respectively, Z is related to the 
efficiency rj of converting the heat current Jq (between the 
baths) into the electric power P which is generated by attach- 
ing a thermoelectric element to an optimal Ohmic impedance. 



Namely, in the linear regime: 



P 



^carnot 



^ZT+l - 1 



VZT +1 + 1 



(2) 



where jycamot = 1 ^ Tc/Th is the Carnot efficiency and 
T = (Th + Tc)/2. Thus a good thermoelectric device is 
characterized by a large value of the non-dimensional figure 
of merit ZT. 

Since the 1960's many materials have been investigated 
but the maximum value found for ZT was achieved for the 
(Bii_2;Sbj;)2(Sei_j,Tey)3 alloy family with ZT « 1. How- 
ever, values ZT > 3 are considered to be essential for thermo- 
electrics to compete in efficiency with mechanical power gen- 
eration and refrigeration at room temperatures. The efforts 
recently focused on a bulk of new advanced thermoelectric 
materials and on low-dimensional materials, and only a small 
increment of the efficiency, ZT < 2.6, has been obtained [3]. 

One of the main reasons for this partial success is a limited 
understanding of the possible microscopic mechanisms lead- 
ing to the increase of ZT, with few exceptions [4]. From a 
dynamical point of view, cross effects in transport have been 
barely studied [5, 6]. So far, the challenge lies in engineering a 
material for which the values of S, a and k can be controlled 
independently. However, the different transport coefficients 
are interdependent, making optimization extremely difficult. 

In this paper we take a completely different approach. In- 
spired by kinetic theory of ergodic gases and chaotic billiards, 
we show that large values of ZT, in principle approaching 
to Carnot's efficiency, can be obtained when the energy of 
the carrier particles does not depend on the thermodynamic 
forces. 

In the linear response regime (see e.g. [7]), one writes 
the general expressions for the heat current Jq and the elec- 
tric cuiTent Je through an homogeneous sample subjected to 
a temperature gradient dxT and a electrochemical potential 
gradient dxll as 



Jq= - k' dxT - TaS d^fJ- , 
Je = — aS dxT — a dxp. ■ 



(3) 
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Here and in what follows, we assume that the transport oc- 
curs along the x-direction and the temperature is given in units 
where the Boltzmann constant fee = 1. 

The electrochemical potential is the sum of a chemical and 
an electric part jl = ^ + He, where /j. is the chemical poten- 
tial of the particles and, if e is the particle's charge, /ig = e<P 
is the work done by the particles against an external electric 
field £ = — 8x4). From (3) the usual phenomenological rela- 
tions follow: if the thermal gradient vanishes, d^T = 0, then 
Je = —a 8x4) = o-£, since for an isothermal homogeneous 
system /j. is uniform. If the electric current vanishes, Je = 0, 
then dxfi = S dxT, which is the definition of the Seebeck 
coefficient, and Jq = —ndxT where k = k! — TaS^ is the 
usual thermal conductivity (see e.g. [9]). 

From the theory of irreversible thermodynamics, n and /ig 
cannot be determined separately; only their combination in /l 
appears in (3) [8J. Based on this equivalence, in what follows 
we take into account the chemical part only, i.e., p. = n. 

Our aim is to study thermoelectricity from an "energy trans- 
port" point of view. To linear order, the energy and parti- 
cle density currents Ju and Jg respectively, can be written in 
terms of the Onsager matrix L [9, 10] as 



Jq — Lqu dx[ — \+ Lgg dx (~^) ) 



(4) 



where Jg = eJg. In the absence of magnetic fields, the On- 
sager reciprocity relations states that L is symmetric, Lug = 
Lgu. From the entropy balance equation for open systems 



Ju — Jq fJ'Jg 



(5) 



and substituting Jq in (3) in favor of J„ and comparing the 
resulting equations with (4) it follows that the transport coef- 
ficients can be written in terms of the i-coefficients as 
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(6) 

Eqs. (4) and (6) are completely equivalent to the description 
(3). Furthermore, from Eq. (1), we obtain for the figure of 
merit 



ZT = 



{Lug ~ t^Lgg) 

dctL 



(7) 



Expressions (4) and (7) provide a very convenient way for 
numerical or analytical evaluation of ZT for different kinds 
of dynamical models. 

The second law of thermodynamics only requires that L is 
positive definite. Therefore, from (7) it is clear that the second 
law does not impose any upper bound on the value of ZT. 
Furthermore, the crucial observation is that the Carnot's limit 
ZT = oo is reached when the energy density current and the 
electric current become proportional, since then det L = 0. 
Suppose for example that both energy and charge are carried 



only by non-interacting particles, like in a dilute gas. Then the 
microscopic instantaneous currents per particle at position x* 
and time t, are 



3u{x\t) = E{t)vx{x{t),t)5{x* - x{t)) , 
je{x*,t) = eVxS{x* - x{t)) , 



(8) 



where E is the energy of the particle, x its position and Vx 
its velocity along the field. The thermodynamic averages of 
the two currents (appearing in Eq. (4)) become proportional 
precisely when the variables E and Vx are un-correlated 

(E) 

Ju = {ju) = {E){Vx} = ^Oe) = {E)Jg . (9) 

Therefore, ZT = oo follows from the fact that the average 
particle's energy {E) does not depend on the thermodynamic 
forces. In the context of classical physics this happens for 
instance in the limit of large number of internal degrees of 
freedom (d.o.f.), provided the dynamics is ergodic. 

We consider an ergodic gas of non-interacting, electrically 
neutral particles of mass m with rfint internal d.o. f. (rotational 
or vibrational), enclosed in a dimensional container. To 
study the non-equilibrium state of such dilute poly-atomic 
gas we consider a chaotic billiard channel (Fig. 1) connected 
through openings of size A to two reservoirs of particles which 
are idealized as infinite chambers with the same poly-atomic 
gas at equilibrium density g and temperature T. From the 
reservoirs, particles are injected into the channel at a rate 7 
obtained by integration over energy of the appropriate canon- 
ical distribution to give 
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(10) 



The particle injection rate 7 is related to the value of the 
chemical potential at the reservoirs which, for polyatomic 
molecules with a total oiD = d + dint d.o.f., reads 



At = Tin 



cp 7 

T(-D+l)/2 



(11) 



where cj^ is a -D-depending constant. Furthermore, averaging 
the energy of the injected particles over the canonical distri- 
bution, denoted as (£"), we obtain the rate at which energy is 
injected from the reservoirs as £ = 7(-E) = ")T{D + l)/2. 

Let pt{l) be the transmission probability of the channel of 
length I. For a billiard system of noninteracting particles [11], 
the density currents J„, Jg assume a simple form: they are 
Pt{l) times the difference between the left and right corre- 
sponding injection rates, £, 7, respectively, namely 

Jg = Pt (7L - 7r) , Ju = Pt (el - £r) . (12) 
Using (1 1) to eUminate 7 in favor of ^ we obtain, 

' (13) 



"^^ ~ (27rm) 
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FIG. 1: The open Lorentz gas system and a typical particle's tra- 
jectory. The composite particle (schematically represented as a 
molecule) is scattered from fixed disks of radius R disposed in 
a triangular lattice at critical horizon, i.e., the width and height 
of the cells are Aa; = 2R and Ay = 2W respectively, where 
W — 47?/-v/3 is the separation between the centers of the disks. The 
channel is coupled at the left and right boundaries to two thermo- 
chemical baths at temperatures Tl and Tr and chemical potentials 
/XL and fi-R, respectively. 



Taking total differentials of (13) in the variables 1/T and fi/T 
and comparing the resulting expression with Eq. (4) we obtain 
exact microscopic expressions for the Onsager coefficients, 
namely 



L 



(27rm)i/2 
Xptil)l 



qT 



3/2 



(14) 



(27rm)i/2 \^ 2 

\pt{l)l {D + l){D + ?,) 
" (27rm)V2 4 ^ " 

Note that for a chaotic billiard channel with a diffusive dy- 
namics, the transmission probability decays as pt{l) oc l^^ 
which means that all the elements of the Onsager matrix L 
become size independent. 

Finally, plugging (14) into (7) and noting that Cy = D/2is 
the dimensionless heat capacity at constant volume of the gas, 
we obtain 



1 



Cy V T 



(15) 



where for simplicity we have called cy = Cy + 1 /2. A partic- 
ular case of (15) was previously obtained, for noninteracting 
monoatomic ideal gases in 3 dimensions [12]. 

In absence of particles' interaction, ZT is independent of 
the sample size I and depends on the temperature only through 
the chemical potential term. This is due to the fact that with no 
interactions, pt depends on the geometry of the billiard only. 
From a physical point of view this means that the mean free 
path of the gas particles is energy independent. Were particles 
interacting, pt would depend on the local density and tem- 
perature of the gas, leading to a more realistic situation [13]. 



Interestingly enough, we have found that Eq. (15) is an up- 
per limit of the interacting case when the interaction strength 
vanishes [14]. 

We shall now confirm Eq. ( 1 5) with a very simple numerical 
demonstration of a 2-dimensional chaotic Lorentz gas channel 
of particles elastically colliding with circular obstacles of ra- 
dius R. In what follows, we fix the unit length setting R = 1. 
The geometry of the model is shown in Fig.l. 

We consider composite particles with dint > 1 internal ro- 
tational d.o.f.. Each "particle" of mass m can be imagined as 
a stack of djnt small identical disks of mass m/dint and ra- 
dius r R, rotating freely and independently at a constant 
angular velocity LUi, i = 1, . . . dint- The center of mass of the 
particle moves with velocity v = {vx,Vy). 

At each collision of the particle with the boundary of the 
billiard (either one of the circular obstacles or the outer wall) 
an energy exchange among all the D d.o.f. occurs according 
to the following collision rules 



1 - Vdhit _ 
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dint(l + '7c'int) 



UJk , 



where {vn,vt) are the normal and tangent components of 
the particle velocity at the collision point and the parameter 
T] = Q/nir^, with Q being the moment of inertia of each 
small fictitious internal disk. The primed (unprimed) quan- 
tities refer to their values after (before) the collision. These 
collision rules are a generalization of the ones introduced in 
[6]. Thus, they are deterministic, time reversible and preserve 
the energy and local angular momentum. The derivation of 
the collision rules will appear elsewhere [14]. 

First we considered a closed system in a finite container 
and we checked energy equipartition among all d.o.f.. Then 
we have opened the system from both ends and allowed it to 
exchange particles with the two baths at temperatures Tl, Jr. 
and with chemical potentials /iL, Mr. The coupling among the 
system and the baths is defined as follows: whenever a parti- 
cle in the system crosses the opening which separates it from 
the bath, it is removed from the system. On the other hand, 
with frequency 7, particles are injected into the system, with a 
velocity distributed according to the canonical distribution at 
the corresponding temperature. 



Pn[Vn) = — |w„|exp I — 



Pt{vt) 



m, 
2^ 



exp 



4L>7rr 
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2T ^ 
mujf \ 
ADTJ 



(17) 



reflecting the assumption that the bath is an ideal gas at equi- 
librium temperature T. 
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FIG. 2: Figure of merit ZT as a function of the number of d.o.f., at 
= 0. For each D, the transport coefficients were obtained from two 
different simulations in a channel of 10 cells, at mean particle density 
Q = 0.446 and mean temperature T = 1000 with: a) /^T /T = 
0.04, A{fi/T) = 0, and 6) AT = 0, A{fi/T) = 0.04. The injection 
rates 7l and 'jr are obtained from A(/i/r), using (11) and (10). The 
dashed line corresponds to {D + l)/2. In the inset, the dependence 
of ZT on the length of the channel / is shown for D = 9. The dashed 
line shows the corresponding expected value. 
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FIG. 3: Probability distribution function of the energy per degree of 
freedom Ed = E/D determined from equilibrium simulation with 
Tl = Tb,, /xl = /xr, for different number of freedoms: D = 3 
(red), D = 13 (green) and D — 43 (blue). The dashed curves are 
the theoretical ("Chi-square" Xd) distributions of Ed- In the inset, 
the corresponding probability distribution functions for the x- com- 
ponent of the velocity P{vx) (blue) and for the angular momentum 
of one of the particle's disks P{ijJi) (green) is shown, for D = 3. 
The dashed curves are the theoretical Gaussian distributions. 



By two different simulations with two linearly indepen- 
dent sets of thermodynamic forces we have numerically de- 
termined the Onsager's matrix L and the value of ZT (see 
Fig. 2). Numerical results excellently confirm the theoretical 
prediction (15). We have also carefully checked that all On- 
sager coefficients, or conductivities, decay with the size I of 
the system as l/l which indicates a diffusive transport. 

The simple mechanism for the growth of ZT with (ii„t is 
nicely illustrated in Fig. 3 which shows that the particle veloc- 
ity Vx has a Maxwellian (Gaussian) distribution (inset), while 
the equilibrium distribution of the particle energy per degree 
of freedom Eo becomes more and more sharply peaked, and 
thus de-correlated from Vx as dint grows. 

In conclusion, we have discovered a simple general theoret- 
ical mechanism which may find a way to applications of ther- 
moelectricity in real world materials. Even though the case of 
an ionized polyatomic gas may seem a little artificial in this 
context, there may be other important instances where each 
charge would be caiTied by many effectively classical d.o.f.. 
We have also performed the first numerical computation of 
ZT from deterministic microscopic equations of motion. Our 
method can easily be implemented for more realistic models 
where also quantum effects can be taken into account. 
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